一 ChinaXiv 合 作 期 刊 


第 51 卷 第 12 期 全 局 学 组 VolL51 No.12 


2015 年 12 月 第 1507-1515 页 ACTA METALLURGICA SINICA Dec. 2015 pp.1507-1515 


耦合 迹 生 的 TW 耻 钢 多 晶体 塑性 变形 行为 研究 
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摘 要 基于 已 建立 的 单 晶 体 塑性 模型 , 建立 了 耦合 迹 生 的 挛 生 诱发 塑性 TWIP) 钢 多 晶体 塑性 模型 ， 人 
晶 均 匀 化 处 理 相 邻 晶 粒 间 的 几何 协调 和 应 力 平衡 条 件 , 获得 了 单 晶体 与 多 晶体 状态 变量 的 关系 , 开发 了 基于 ABAQUS/ 
UMAT 的 计算 程序 . 采用 EBSD 研 究 了 TWIP 钢 拉 伸 应 变 分 别 为 0.27 和 0.6 时 的 织 构 变 化 , 并 对 模型 进行 了 应 力 应 变 及 织 构 
演化 的 验证 . 用 该 本 构 模 型 分 别 建立 了 拉 伸 、 压 缩 和 扭转 3 种 简单 加 载 条 件 下 的 有 限 元 模型 , 分 析 了 不 同 变形 条 件 下 的 宏观 
力学 响应 及 织 构 演化 规律 . 结果 表明 : 拉 伸 变形 过 程 中 , 应 变 硬 化 现象 和 织 构 密度 水 平 随 应 变 增加 而 增强 ; 在 压缩 过 程 中 ， 

织 构 类 型 随 应 变 增加 而 发 生变 化 , 但 是 织 构 密度 水 平 基本 不 变 ; 而 在 扭转 过 程 中 , 当 扭 转 应 变 较 小 时 , 基本 无 织 构 形成 , 随 
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本 一 着 应 变 增加 , 织 构 逐 渐 显 现 出 来 , 这 是 因为 变形 较 小 时 , 圆柱 沿 径 向 方向 内 部 变形 量 较 小 , 故 织 构 不 明显 . 
关键 词 TWIP 钢 , 晶体 塑性 , 多 晶 均匀 化 , 织 构 预测 , 塑性 变形 
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ABSTRACT Twinning induced plasticity (TWIP) steel exhibits high strength and exceptional plasticity due to 
the formation of extensive twin under mechanical load and its ultimate tensile strength and elongation to failure 
ductility-value can be as high as 5x10’ MPa.%, which provide a new choice for automobile in developing the light- 
weight and improving safety. Generally, due to the texture was formed during process of plastic deformation, metal 
material appear anisotropic behavior. The deformation mechanisms, responsible for this high strain hardening, are 
related to the strain-induced microstructural changes, which was dominated by slip and twinning. Different defor- 
mation mechanisms, which can be activated at different stages of deformation, will strongly influence the stress 


strain response and the evolution of the microstructure. In this work, to predict the texture evolution under different 
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loading conditions and understand these two deformation mechanisms of plastic deformation process, a polycrystal 


plasticity constitutive model of TWIP steel coupling slip and twinning was developed based on the crystal plastici- 


ty theory and single crystal plasticity constitutive model. A polycrystal homogenization method to keep geometry 


coordination and stress balance adjacent grains was used, which connected the state variables of single crystal and 


polycrystal. And then the model was implemented and programed based on the ABAQUS/UMAT platform. The tex- 


ture evolution was obtained by EBSD at strain 0.27 and 0.60, respectively. The finite element models of tensile, 


compression and torsion processes were built by using the constitutive model. The mechanical response and texture 


evolution during plastic deformation process of TWIP steel were analyzed. The results show that with the increas- 


ing of the strain, the Strain hardening phenomenon and texture density enhanced during the tensile process. Al- 


though texture types changed, texture density unchanged during the compression process. Owing to deformation in- 


creasing along the diameter direction, there is no obvious texture inside the cylinder when torsion deformation is 


small, texture emerged and enhanced gradually with the increasing of strain. 


KEY WORDS TWIP steel, crystal plasticity, polycrystal homogenization, texture prediction, plasticity deforma- 


tion 


挛 生 诱发 塑性 (twinning induced plasticity, 
TWIP) 钢 不 仅 具 有 高 抗 拉 强 度 和 高 硬化 率 , 而 且 其 
有 优良 的 塑性 、 万 性 和 成 形 性 能 ,为 现代 汽车 发 展 
轻 量化 和 提高 安全 性 等 方面 提供 了 新 的 选择 , 引起 
国内 外 研究 者 的 广泛 关注 号. 一 般 来 说 , 经 过 变形 后 
的 多 晶体 金属 材料 通常 出 现 各 向 异性 行为 , 这 主要 
由 于 多 晶体 在 塑性 变形 过 程 中 形成 变形 织 构 中 . 
TWIP 钢 塑性 变形 是 由 位 错 滑 移 和 杰 生 变形 2 种 机 
制 共 同 作用 的 结果 , 当 唱 体内 挛 晶 体积 达到 其 饱和 
值 时 , 引起 唱 格 发 生 转 向 , 转向 后 位 错 滑 移 在 新 的 
滑 移 面 上 继续 进行 . 因此 , TWIP 钢 织 构 的 形成 是 位 
昔 滑 移 和 挛 生 共同 作用 产生 的 中 . 材料 变形 织 构 的 
形成 与 加 载 条 件 相关 , 通过 对 不 同 变 形 条 件 下 织 构 
的 预测 , 控制 材料 的 成 形 工艺 以 实现 材料 成 形 性 能 
最 优化 . 
研究 人 员 ® 采用 现代 显 微 测试 技术 和 理论 建 
模 对 TWIP 钢 进行 了 研究 , 解释 了 其 塑性 变形 是 由 
挛 生 和 滑 移 共 同 作 用 引起 的 ; 并 进一步 分 析 了 滑 移 
和 挛 生 的 作用 规律 , 为 定性 分 析 滑 移 和 挛 生 对 宏观 
塑性 变形 的 贡献 提供 了 参考 , 但 没有 给 出 各 挛 生 系 
和 滑 移 系 的 演化 规律 . 本 文 作者 前 期 工作 "对 单 晶 
体 塑 性 模型 研究 得 到 了 TWIP 钢 单 晶 体 塑性 变形 过 
程 中 各 挛 生 系 ` 滑 移 系 的 演化 过 程 及 挛 生 对 位 错 滑 
移 的 影响 , 但 不 能 用 于 描述 多 晶体 织 构 的 演化 . 此 
外 , 研究 人 员 以 单 唱 体 塑 性 模型 为 基础 , 建立 了 
Sachs 模型 吓 、Taylor 模 型 号、 自治 (self consistent SC) 
模型 中 和 有 限 元 模型 "等 多 晶体 本 构 模 型 , 对 多 唱 
体 进行 了 变形 织 构 的 实验 和 模拟 研究 . 目前 , 由 于 
Sachs 模型 假设 多 晶 变 形 时 各 品 粒 所 承受 的 应 力 状态 
与 宏观 张 量 一致 , 因此 在 各 唱 粒 唱 界 之 间 产 生 的 应 变 


不 连续 , 与 实验 结果 不 太 吻 合 ; Taylor 模 型 满足 了 唱 界 
侧 的 应 变 协调 条 件 , 与 实验 结果 吻合 较 好 ,但 在 
晶 界 处 的 应 力 平衡 却 不 能 得 到 满足 ; 自 洽 模型 在 应 
复杂 本 构 关 系 时 , 需要 反复 近代 计算 , 应 用 范 
也 比较 有 限 . 之 后 , 研究 人 员 "”" 将 晶体 塑性 理论 与 
大 型 商用 有 限 元 软件 结合 , 在 介 观 尺度 上 形成 了 一 
种 新 的 织 构 模 拟 方法 , 即 晶 体 塑性 有 限 元 法 (crystal 
plasticity finite element method, CPFEM), 该 方法 不 
仅 能 够 模拟 变形 过 程 中 织 构 的 演化 规律 , 而 且 能 

于 与 织 构 相 关 的 各 种 性 能 响应 . 同时 , 结合 电子 
背 散 射 衍 射 (electron back-scattered diffraction, EB- 
SD) 分 析 技 术 , 能 够 实时 、 快 速 和 统计 性 地 获取 多 唱 
体 中 各 个 晶 粒 的 取向 信息 . 研究 人 员 W”>3 采 用 晶体 
塑性 理论 , 结合 有 限 元 模拟 不 同 变 形 条 件 下 的 织 构 
演化 , 得 到 了 与 实验 一 致 的 良好 结果 , 但 并 未 对 李 
生 和 滑 移 对 织 构 演 化 的 影响 做 出 定量 分 析 ; Li 等” 
对 TWIP 钢 滑 移 和 杰 生 对 应 变 硬化 和 织 构 演 化 的 贡 
献 做 了 定量 分 析 , 认为 在 初始 变形 阶段 主要 依靠 这 
生变 形 , 随 着 应 变 的 增加 , 挛 生 对 塑性 变形 的 贡献 
很 快 被 滑 移 取 代 , 最 终 挛 生 和 滑 移 的 贡献 达到 平 
衡 , 但 并 未 对 不 同 加 载 条 件 对 织 构 形成 的 影响 进行 
深入 分 析 . 

本 工作 对 退火 处 理 后 的 Fe-22Mn-0.6C 钢 (质量 
分 数 , %) 进 行 单 轴 加 载 实验 , 获得 拉 伸 条 件 下 宏观 
力学 响应 及 织 构 演化 特点 ; 基于 已 建立 的 TWIP 钢 
单 晶 体 塑性 本 构 模 型 % 采用 多 唱 均 匀 化 处 理 建 立 
联系 单 唱 体 与 多 晶体 状态 变量 的 多 晶体 塑性 本 构 
模型 , 对 该 模型 在 ABAQUS/UMAT 平 台 上 进行 二 
次 开发 并 进行 数值 模拟 , 结合 实验 确定 模型 中 相关 
参数 并 验证 模型 可 靠 性 ; 进一步 将 上 述 本 构 模 型 应 
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用 于 其 他 简单 加 载 条 件 (压缩 和 扭转 ) 下 的 应 力 应 变 
关系 及 织 构 演 化 , 分 析 不 同 加 载 条 件 对 应 变 硬 化 的 
影响 及 织 构 形成 的 贡献 . 
1 实验 方法 

本 工作 实验 用 TWIP 钢 多 晶体 材料 的 主要 化 学 
成 分 为 Fe-22Mn-0.6C. 原材料 采用 Ar 气氛 保护 , 电磁 
感应 炉 真 空 熔炼 , 浇铸 成 板 坏 , 板 坏 经 热 轧 至 3.5 mm 
然后 冷 轧 至 1.5 mm, 在 SRJX-8-13 箱 式 电阻 炉 中 进 
行 退 火 处 理 , 温度 为 800“C, 保温 时 间 75 min. 在 二 
辊 冷 轧 试验 机 上 经 过 多 道 次 轧 制 后 厚度 为 1 mm, 再 
在 四 辊 冷 轧 试验 机 上 经 过 多 道 次 轧 制 最 终 获 得 厚 
度 为 0.5 mm 的 TWIP 钢 钢板 . 为 了 消除 材料 的 初始 
轧 制 织 构 , 对 TWIP 钢 钢板 进行 均匀 化 退火 处 理 , 温 
度 为 800 'C, 保温 时 间 240 s™. 

为 了 掌握 TW 了 PP 钢 多 晶体 在 不 同 变形 条 件 下 的 
应 力 应 变 特点 及 织 构 演化 规律 , 对 经 退火 处 理 后 的 
TW 了 PP 钢 钢板 沿 轧 制 方向 取样 , 试 样 在 MTS810 材料 
试验 机 上 进行 2 组 不 同 应 变量 的 拉 伸 实验 , 拉 伸 速 
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数 ， 广 为 第 5 挛 生 系 的 挛 晶 体积 分 数 ， 
上 述 模 型 描述 了 TWIP 钢 单 晶 体 塑 性 变形 过 程 
中 各 挛 生 和 滑 移 系 的 应 变 硬化 及 其 演化 特点 , 而 实 
际 金属 材料 多 为 多 晶体 结构 , 因而 需要 通过 多 唱 均 
匀 化 处 理 , 建立 单 晶 体 状 态 变量 与 多 晶体 状态 变量 
间 的 关系 , 使 2 个 相 邻 的 晶体 在 物理 条 件 上 满足 几 
可 协调 和 应 力 平衡 条 件 . 
Roters 等 中 基于 连续 介质 理论 提出 了 有 限 元 
均匀 化 模型 , 假设 一 个 单元 代表 一 个 晶体 , 材料 质 
点 的 应 力 响应 由 单个 晶体 的 应 力 响应 决定 . 为 了 
将 文献 [9] 中 的 单 晶 变量 参数 与 多 蝇 变 量 参 数 进 
行 区 别 , 在 未 进行 特殊 说 明 时 , 将 多 晶 变 量 参数 上 
方 加 “一 ” 
根据 连续 介质 大 变形 理论 , 未 变形 时 多 晶体 质 
点 的 参考 构 形 B 由 材料 的 初始 状态 决定 , 施加 外 载 
答 变形 后 多 晶体 的 宏观 塑性 变形 梯度 中 为 : 
F=V,y (4) 


PE 


率 均 为 8 mm/min. 第 一 组 拉 伸 后 测 得 真 应 变 为 
0.27; 第 二 组 拉 伸 至 试 样 断裂 , 断后 测 得 真 应 变 为 
0.60. 将 未 拉 伸 试 样 和 拉 伸 后 的 2 组 试 样 再 次 取样 
进行 金 相 研磨 并 电解 抛光 , 采用 EBSD 技术 观察 织 
构 状 态 . 
2 多 晶体 塑性 均匀 化 方法 和 数值 模拟 
2.1 多 晶体 塑性 的 均匀 化 方法 

文献 [9] 中 给 出 了 单 晶体 塑性 本 构 关 系 , 分 析 了 
单 晶 体 塑 性 变形 过 程 中 滑 移 系 和 挛 生 系 的 演化 规 
律 , 并 引入 村 晶体 积分 数 饱 和 值 , 揭示 了 单 晶体 内 
挛 晶 体积 的 演化 规律 及 挛 生 对 位 错 滑 移 的 影响 , 滑 
移 和 挛 生 的 硬化 模 量 可 用 如 下 形式 表征 入 

$i a 
Sra = Sa tl 2 了 (1) 


式 中 , s,,,。 和 ,分 别 为 第 n+tl1 步 和 第 n 步 计算 时 滑 
移 阻力 , n 为 计算 步 数 ; ye 为 第 a 滑 移 系 的 滑 移 速 
率 ; At 为 时 间 增 量 ; h., 和 s, 分 别 为 第 a 滑 移 系 的 
硬化 率 和 滑 移 阻 力 饱 和 值 , 可 以 表示 为 : 


认 中 + by 中 (2) 
B 


Ss SD bg (3) 
B 


式 中 , h, 和 s。 分 别 为 未 发 生 挛 生 时 的 清 移 系 的 硬 
化 率 和 滑 移 阻力 饱和 值 ，s, , c, 5b 都 是 材料 硬化 参 


式 中 , 不 为 宏观 塑性 变形 梯度 , 表示 多 晶体 弹性 和 
塑性 变形 梯度 的 总 和 ; 达 为 多 晶体 材料 质点 矢 径 ， 
表示 从 坐标 原点 到 物质 点 的 矢量 ; 了 为 多 晶体 材料 
质点 的 非 线性 变形 量 ，Y.7 为 了 对 元 求 偏 导 . 

严 可 以 分 解 为 Ps: 


F=F.:F, (5) 
式 中 , 扩 为 多 晶体 弹性 变形 ( 唱 格 畸变 ) 和 刚性 转动 
的 合成 变形 梯度 , 屎 ,为 多 晶体 滑 移 和 挛 生 变形 引起 
的 塑性 变形 梯度 . 

多 晶体 弹性 变形 (品格 畸变 ) 和 刚性 转动 的 合 
成 变形 梯度 、 滑 移 和 这 生变 形 引 起 的 塑性 变形 梯度 
均 可 表示 为 与 所 有 唱 体 相对 应 变形 梯度 的 体积 平 
均值 : 


m1 

=| (6) 

| 

4 0) 
v=Jdv (8) 


式 中 , FF. 为 单 晶 体 弹 性 变形 ( 晶 格 畸变 ) 和 刚性 转动 
的 合成 变形 梯度 ,，F, 为 单 晶 体 滑 移 和 杰 生 变形 引起 
的 塑性 变形 梯度 , 8 为 变形 后 多 晶体 的 构 形 , Y 为 多 
晶体 体积 . 
由 文献 [9] 可 知 , 多 晶体 第 二 Piola-Kirchhoff 平 

均 应 力 与 到 的 关系 为 : 
T=F'. {det(F)T}: F." (9) 
式 中 , TT 为 多 晶体 第 二 Piola-Kirchhoff 平 均 应 力 , 表 
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示 所 有 晶体 第 二 Piola-Kirchhoff 应 力 的 体积 平均 
值 ; 无 ' 为 素 的 道 矩 阵 ; det( 瓦 ) 为 到 的 行列 式 ; 7 
为 多 晶体 Cauchy 平均 应 力 , 表示 所 有 晶体 Cauchy 
应 力 的 体积 平均 值 ; 严 " 为 到 转 置 后 的 逆 和 矩阵 . 

而 多 晶体 第 一 Piola-Kirchhoff 平 均 应 力 与 开关 
系 为 : 


五 = 到 .不 (10) 
五 = $lPav (11) 


式 中 , P 为 多 晶体 第 一 Piola-Kirchhoff 平 均 应 力 , 表 
示 所 有 晶体 第 一 Piola-Kirchhoff 应 力 的 体积 平均 
值 ; P 为 单 晶体 第 一 Piola-Kirchhoff 应 力 . 
上 述 均匀 化 模型 中 , 每 个 有 限 单元 代表 一 个 唱 
体 , 单元 的 边界 正好 是 晶体 的 边界 . 而 有 限 元 中 绝 
大 多 数 单元 是 协调 元 , 在 单元 的 边界 上 满足 几何 协 
调 条 件 , 同时 也 在 弱 形 式 上 满足 应 力 平 衡 条 件 . 那 
么 相 邻 晶体 之 间 的 几何 协调 和 应 力 平衡 条 件 就 可 
以 得 到 满足 . 
2.2 多 晶体 塑性 的 数值 模拟 过 程 

本 工作 中 多 晶体 塑性 的 数值 模拟 过 程 如 图 1 所 
示 . 利用 ABAQUS/UMAT 将 初始 时 刻 所 有 晶体 进 
行 随机 取向 分 布 , 并 定义 所 有 唱 体 的 材料 参数 和 计 
算 参 数 ; 利用 单 晶体 塑性 本 构 模型 计算 每 个 晶体 在 
变形 过 程 中 的 状态 变量 , 同时 根据 多 晶 均 匀 化 处 理 
方法 调节 相 邻 晶体 间 应 力 平衡 条 件 和 几何 协调 关 
系 ; 计算 完成 后 输出 所 需 状 态 变 量 . 
3 耦合 迹 生 的 多 晶体 塑性 模型 参数 确定 及 实 

验 验 证 

3.1 多 晶体 塑性 模型 的 参数 确定 
结合 文献 [9,10] 中 TWIP 钢 单 晶体 塑性 本 构 模 


1 多 晶体 塑性 本 构 模型 数值 模拟 的 流程 图 


型 中 的 参数 获取 方法 , 并 取 初 始 试探 值 经 过 不 断 修正 
可 得 到 本 构 模 型 中 的 材料 参数 . 四 阶 弹性 张 量 C 中 的 
3 个 材料 常数 分 别 取 为 Cu=198 GPa, Cw=125 GPa, 
Cu=122 GPa; 挛 生 体积 分 数 阔 值 有 取 为 0.2; 在 准 静 
态 变形 时 参考 塑性 剪 切 率 y, 取 为 0.001 s'; 率 敏 感 
系数 m=0.0074; 计生 和 滑 移 阻 力 比值 6=1.16, 其 余 
参数 见 文 献 [9]. 
3.2 多 晶体 拉 伸 过 程 的 有 限 元 模型 

本 工作 通过 ABAQUS/UMAT 二 次 开发 对 
TWIP 钢 多 晶体 塑性 本 构 模 型 进行 数值 模拟 , 并 进 
一 步 建立 TWIP 钢 多 晶体 拉 伸 过 程 的 有 限 元 模型 ， 
以 再 现 TWIP 钢 多 蝇 试 样 单 向 拉 伸 条 件 下 宏观 应 力 
应 变 关系 及 织 构 演化 状态 , 如 图 2 所 示 . 该 模型 考虑 
实际 拉 伸 过 程 中 试 样 有 效 长 度 范围 内 加 载 、 变 形 和 
初始 织 构 分 布 特点 , 可 在 其 有 效 长 度 范 围 内 取 一 特 
征 段 进行 建 模 , 如 图 2a 所 示 的 红色 区 域 , 特征 区 域 
尺寸 为 2.5 mmx0.25 mmxl mm: 考虑 到 特征 段 的 对 
称 特点 , 左 端 和 后 面 分 别 关 于 xz 面 和 zy 面 对 称 约 
束 , 下 表面 z 方 向 约束 , 上 表面 给 以 z 方 向 正 向 位 移 ， 
为 保证 计算 的 收敛 性 , 网 格 单元 采用 C3D8 单元 (三 
维 8 节点 实体 单元 ), 均匀 划分 单元 数量 为 320, 如 图 
2b 所 示 ; 所 得 某 一 晶体 初始 取向 信息 如 图 2c 所 示 . 
图 3a 为 TWIP 钢 多 上 晶体 经 退火 处 理 后 未 变形 时 
的 极 图 . 可 以 看 出 , 经 过 退火 后 的 TWIP 钢 在 {100}， 
{110} 和 {111} 面 上 没有 出 现 明显 的 退火 织 构 , 3 个 极 
图 中 的 极 密度 比较 均匀 且 较 低 , 晶 粒 取向 为 任意 分 
布 . Barbier 等 所 所 研究 的 TWIP 钢 经 过 退火 处 理 后 
仍 有 明显 织 构 , 因此 考虑 了 初始 织 构 对 TWIP 钢 拉 
伸 变 形 过 程 中 织 构 演 化 的 影响 . 本 工作 忽略 初始 退 
火 织 构 对 TWIP 钢 多 唱 材 料 影 响 , 认为 经 过 退火 处 


State variables 
Deformation gradient F,, 
Inverse matrix of plastic deformation gradient F 
Secondary PK stress 也 
Slip resistance s, 


pry 


Deformation gradient F, 
Schmid tensor su and s, 
Material constant 
Step length 


Fig.1 Schematic representation of polycrystal constitutive model 
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(a) (b) (c) 
Color online 


图 2 挛 生 诱发 塑性 TWIP) 钢 多 晶体 拉 伸 过 程 的 有 限 元 模型 
Fig.2 Finite element model of macroscopic tensile process for twinning induced plasticity (TWIP) steel by using CPFE 


model (unit: mm) 


(a) schematic of macroscopic tensile specimen 
(b) finite element model of polycrystal 
| ) 


(c) initial orientation of one crystal after rotation certain Euler angle 
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Color online 


图 3 退火 处 理 和 晶 粒 取向 随机 分 布 后 的 TWIP 钢 在 {100}, {110} 和 {111} 面 织 构 
Fig.3 Experimental textures after annealing treatments (a) and simulated textures after grains random orientation (b) in 


{100}, {110} and {111} planes for TWIP steel (RD 一 rolling direction; TD 一 transverse direction) 
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理 后 所 有 唱 粒 取向 均匀 分 布 . 结合 初始 条 件 下 的 实 。 生 新 的 {111}<110> 织 构 , 其 密度 水 平均 在 8 左右 . 同 
验 结果 , 对 所 有 晶 粒 进行 随机 取向 分 布 , 所 有 晶 粒 。 时 随 着 拉 伸 应 变 的 增加 , 拉 伸 织 构 的 锋 锐 程度 增 
随机 取向 分 布 后 织 构 如 图 3b 所 示 . 可 以 看 出 , 晶 粒 加. 通过 对 比 结果 可 知 , 该 模型 能 较 好 地 预测 织 构 
随机 取向 分 布 后 的 TWIP 钢 在 {100}, {110} 和 {111} 的 演化 , 具有 一 定 可 靠 性 , 可 在 实际 的 TWIP 钢 多 晶 
面 上 没有 出 现 明显 的 织 构 , 晶 粒 取向 为 任意 分 布 ， 材料 塑性 加 工 过 程 中 进行 应 用 . 
与 退火 条 件 下 的 晶体 织 构 有 良好 的 一 致 性 . 4 不 同 加 载 条 件 下 TWIP 钢 多 晶体 塑性 变形 
3.3 宏观 应 力 应 变 和 晶体 取向 的 验证 行为 

本 工作 结合 TWIP 钢 多 品 体 拉 伸 实验 和 所 建立 为 了 进一步 分 析 不 同 加 载 条 件 下 TWIP 钢 多 星 


的 拉 伸 过 程 有 限 元 模型 , 获得 拉 伸 实验 的 应 力 .应 ed 
前 性 变形 特点 及 织 构 的 形成 和 演化 规律 ,基于 所 建 
变 硬化 率 和 这 晶体 积分 数 随 应 变 的 演化 与 其 模拟 体 理 性 到 形 翌 所 及 织 得 的 形成 和 演化 规律 , 基于 所 建 


立 的 多 晶体 塑性 元 模型 , 分 别 瓜 
有 如 辣 4 所 未 可 以 看 由 在 应 灾 为 027 和 00。 人 直人 
、 日 日 Dts | 和 。 一 J 
条 件 下 模拟 获得 的 应 力 应 变 曲 线 与 实验 获得 的 结 


有 限 元 模型 相似 , 所 有 网 格 单元 采用 C3D8 单 元 ， 
压缩 变形 过 程 数值 建 模 采 用 与 上 述 拉 伸 变 形 模型 相 
同 的 特征 段 , 尺寸 均 为 2.5 mmx0.25 mmxl mm, 均匀 
划分 为 320 个 有 限 单 元 , 位 移 边界 条 件 与 拉 伸 变 形 
a 数值 模型 相 同 ， 位 移 加 载 与 拉 伸 变形 方向 相反 , 通 
i a 1 过 加 载 不 同位 移 量 来 控制 变形 量 的 大 小 . 扭转 变形 

ee ee ) 由 机 网 有 限 元 模型 采用 直径 1 mmxl mm 圆柱 体 模 型 , 自由 


同 定 , 在 另 一 端面 施加 扭转 变形 , 通过 施加 不 同 捏 


由 


果 有 较 好 的 一 致 性 ; TWIP 钢 有 较 高 的 强度 和 较 好 
的 塑性 , 抗 拉 强度 达到 1200 MPa, 断后 延伸 率 达 到 
80%, 随 着 应 变 的 不 断 增加 , 应 变 硬化 增强 , 应 变 硬 
化 率 逐 渐 降 低 ; 研究 人 员 芋 % 采用 TEM 观察 了 


型 


AS 


图 5 为 真 应 变 为 0.27 和 0.6 时 的 织 构 图 和 取向 
分 布 函数 (orientation distribution function, ODF) 图 的 


转角 度 来 控制 变形 量 的 大 
图 6 给 出 了 压缩 和 扭 


小 . 
转变 形 过 程 中 应 力 、 硬 化 
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实验 和 模拟 结果 . 可 以 看 出 , 随 应 变 的 增加 TWIP 钢 
多 晶体 逐渐 形成 织 构 , 且 织 构 强 度 不 断 增强 . 从 
ODF 图 实验 和 模拟 结果 可 以 看 出 , 单 向 拉 伸 后 的 
TWIP 钢 多 晶体 主要 是 产生 黄 铜 R 型 织 构 , 包括 


率 和 挛 品 体积 分 数 随 应 变 演化 的 模拟 结果 . 可 以 看 
出 , 压缩 和 扭转 过 程 中 随 应 变 的 增加 应 力 和 这 晶体 
职 分 数 逐 渐 增 加 , 应 变 硬化 率 逐 渐 降 低 . 对 比 图 4 和 
6 可 以 看 出 , 在 拉 伸 和 压缩 变形 过 程 中 , 在 相同 变形 


{111}<112> 织 构 和 {111}<110> 织 构 . 其 中 , 应 变 为 量 条件 下 , 由 于 材料 的 拉 伸 应 变 硬化 模 量 小 于 压缩 
0.27 时 , 主要 为 {111}<112> 织 构 , 其 密度 水 平 在 4 左 。 应变 硬化 模 量 , 致使 压缩 应 力 始 终 大 于 拉 伸 应 力 . 
右 . 当 应 变 为 0.6 时 , {111}<112> 织 构 增 强 , 同时 产 图 7a~d 给 出 了 沿 轧 制 方向 压缩 应 变 为 0.8 和 
1500 . A Exp.resultof stress, £0.27 0.24 -6000 
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4 实验 和 模拟 获得 的 应 力 \ 应 变 硬化 率 及 挛 唱 体积 分 数 的 演化 对 比 


Fig.4 Simulated and experimental results of stress, strain hardening rate and volume fraction of twin (Fr—tension) 
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图 5S 真 应 变 为 0.27 和 0.6 时 的 织 构图 和 ODF 图 的 实验 和 模拟 结果 


Fig.S Experimental (a, e) and simulated (c, g) textures for TWIP steel and corresponding ODFs (b, d, f, h) while strain 6= 
0.27 (a~d) and e=0.6 (e~h) along rolling direction 
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图 6 压缩 和 扭转 过 程 应 力 、 硬 化 率 和 杰 晶 体积 分 数 随 应 变 的 演化 


Fig.6 Evolution curves of stress, strain hardening rate and volume fraction of twin during compression and torsion pro- 


cesses (Fc—compression; 7 一 torsion) 
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图 7 TWIP 钢 压缩 应 变 分 别 为 0.8 和 1.52 以 及 扭转 应 变 分 别 为 0.72 和 1.57 时 模拟 的 极 图 和 ODF 图 
Fig.7 Simulated compression with e=0.8 (a) and e=1.52 (c), and torsion with e=0.72 (e) and 6=1.57 (g) textures for TWIP 


steel along rolling direction, respectively, and corresponding ODFs (b, d, f, h) 


1.52 时 的 {100}, {110} 和 {111} 面 上 的 极 图 和 ODF 
图 . 可 以 看 出 , 单 向 压缩 后 的 TWIP 钢 主要 产生 黄 钢 
和 A 织 构 , 包括 {110}<112> 织 构 和 {110} <111> 织 
构 . 其 中 , 应 变 为 0.8 时 , 主要 为 {110}<112> 织 构 , 其 
密度 水 平 在 3 左右 . 当 应 变 为 1.52 时, {110}<111> 织 
构 增 强 , 其 密度 水 平均 在 3 左右 , 随 着 应 变 增加 , 虽 
然 织 构 类 型 有 变化 , 但 是 织 构 密度 水 平均 较 低 . 
7e~h 给 出 了 绕 轧 制 方向 所 在 轴 扭 转 后 的 
{100}, {110} 和 {111} 面 上 的 极 图 和 ODF 图 . 可 以 看 
出 , 在 扭转 量 较 小 时 , TWIP 钢 多 晶体 几乎 没有 织 构 
形成 , 主要 因为 圆柱 体 扭转 变形 时 材料 变形 量 沿 径 
向 方向 逐渐 增加 , 变形 较 小 时 , 材料 内 部 几乎 没有 
变形 , 多 晶体 内 部 接近 初始 状态 , 故 织 构 状态 不 明 
显 ; 随 着 变形 量 的 增加 , 多 唱 材料 内 部 变形 增加 , 织 
构 逐 渐 明 显 . 而 在 拉 伸 和 压缩 变形 条 件 下 , 所 有 唱 


体 都 将 受到 外 力 而 进行 变形 , 故 即使 在 小 变形 条 件 
下 , 多 晶 材 料 形 成 织 构 也 较 明 显 . 
5 结论 

(1) 成 功 预测 了 拉 伸 过 程 中 黄 铜 R 型 织 构 的 生 
成 , 同时 发 现 随 着 拉 伸 应 变量 的 增加 , 拉 伸 织 构 的 
锋 锐 程度 增加 , 与 实验 结果 有 较 好 的 一 致 性 , 验证 
了 模型 的 可 靠 性 . 

(2) 挛 生 诱 发 塑性 钢 多 唱 材 料 的 拉 伸 应 变 硬化 
模 量 小 于 压缩 应 变 硬化 模 量 , 所 以 在 相同 变形 条 件 
下 材料 的 压缩 应 力 始 终 大 于 拉 伸 应 力 ; 同时 , 拉 伸 
过 程 中 , 随 着 应 变 增 加 , 主要 织 构 类 型 发 生变 化 , 平 
均 密度 水 平 增加 ; 压缩 过 程 中 , 随 应 变 的 增加 , 虽然 
织 构 类 型 变化 , 但 是 平均 密度 水 平 基本 不 变 . 

(3) 圆柱 体 扭转 变形 时 材料 变形 量 沿 径 向 方向 
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逐渐 增加 , 当 扭 转变 形 较 小 时 , 由 于 材料 内 部 变形 较 
小 , 织 构 形成 不 明显 , 随 着 应 变 增 加 , 织 构 逐 渐 明 显 . 
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